Evidence for Landau's critical velocity in superfluid helium nanodroplets from wave 

packet dynamics of attached potassium dimers 
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Femtosecond pump- probe spectroscopy has been used to study vibrational dynamics of potassium 
dimers attached to superfluid helium nanodroplets. Comparing the measured data with theoretical 
results based on dissipative quantum dynamics we propose that the most important effect of the 
helium environment is a general damping of the vibrational dynamics as a result of the interaction 
between dimer and collective degrees of freedom of the helium droplet. The calculations allow us to 
explain crucial experimental findings that are unobserved in gas-phase measurements. Remarkably, 
best agreement with experiment is found for a model where we neglect damping once a wave packet 
moves below a critical velocity. In this way the results provide first direct evidence for the Landau 
critical velocity in superfluid nanodroplets. 

PACS numbers: 33.20.Tp, 67.25.dj, 31.70.Hq, 03.65. Yz, 82.20. Wt, 36.40.-c 



Superfluidity in finite-sized quantum systems is a fun- 
damental issue of current interest, both in liquid helium 
(He) 0, S i and in ultracold atomic gases [3|. As 
first discussed by Landau, superfluidity manifests itself 
through the fascinating effect of frictionless flow below a 
critical velocity In this Letter we present direct evi- 
dence for the existence of a critical velocity in superfluid 
He nanodroplets. 

Spectroscopy of probe molecules embedded inside He 
clusters and He nanodroplets (10 3 — 10 4 atoms) has pro- 
vided some insight into their superfluid behaviour by 
studying the response to external perturbations [3, 0, Q • 
In contrast to pure spectroscopic studies reported ear- 
lier, we present femtosecond pump-probe measurements. 
This technique is well-established for studying vibra- 
tional wave packet (WP) dynamics of diatomic systems 
d, 0, B S EH • Here, we investigate potassium molecules 
(K2) attached to He nanodroplets (He nanodroplet isola- 
tion (HENDI) spectroscopy) He droplets provide a 
versatile test bed for studying relaxation dynamics of the 
immersed species which are cooled to the droplet temper- 
ature (o.38 k) mm. 

The observed real time pump-probe spectra of K2 on 
He droplets differ significantly from previously obtained 
gas phase results [TJ 0|. In this paper we argue that 
dissipative quantum dynamics, here employed through a 
quantum master equation, is crucial for the understand- 
ing of these measurements. Moreover, the dissipation 
has to be combined with stochastic desorption of the 
molecule and shifts of the potential energy surfaces in- 
duced by the He droplets to correctly describe the experi- 
mental results. Remarkably, we find best agreement with 
the experimental findings when damping is omitted for 
very slowly moving WPs. This demonstrates the poten- 
tial of real time studies of vibrational motions to investi- 
gate the Landau critical velocity on the microscopic scale, 
very much in the spirit of vibrational wire resonators [IH 



or quartz tuning forks [16j that probe on larger scales. 

In the pump-probe spectroscopy of K2, potassium 
atoms are attached to He nanodroplets (~ 5 000 atoms) 
which are created in a supersonic expansion of He gas at 
cryogenic conditions Doping conditions are chosen 
such that on average two potassium atoms stick to one 
droplet forming K2 dimers which are weakly bound on 
the surface of the droplets. In a one color pump-probe ex- 
citation scheme the resulting photo ions are recorded 
mass selectively. A laser having 110 fs pulse width has 
been used providing 16 nJ pulses mildly focused in the 
interaction region with the molecular beam. The (3)- 
photon excitation schemes which are most relevant are 
depicted in Fig. [TJ 

The time dependent Schrodinger equation for the state 
ijr(t) = (ip x (t), tpA(t), ■ ■ • ) is solved f ully numerically, as 
detailed for the gas phase in Refs. d,[I3- Here, ipi(t) de- 
notes the vibrational WP in electronic state i. Transition 
dipole moments are assumed to be constant (Condon ap- 
proximation) . In a first step we include He- induced shifts 
of potential energy surfaces (PES). Systematic numerical 
tests reveal that best agreement with experiment is ob- 
tained when only shifting the 2H surface (by — 50cm _1 ) 



which agrees well with findings in [111 . 1 1 81 ] . As the excess 
energy of photo electrons is not observed here, a shift 
of the ionic surface as reported in [l9[ has no effect on 
our results. Next, we take into account dimer desorption 
from the droplets. We choose a constant desorption rate 
leading to a probability density P(t') ~ exp{— t'/ru} for 
the desorption time t' after electronic excitation. Upon 
desorption, the He-induced shifts are set to zero. We find 
that a description solely based on energy shift (parame- 
ter An) and desorption (parameter td) cannot reproduce 
crucial features of the HENDI measurements [201 ] . 

The new essential ingredient of our description is the 
inclusion of damping of vibrational WPs. Dissipation 
of vibrational energy occurs due to the interaction with 
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FIG. 1: Potential curves of the K2 molecule and excitation 
schemes at A = 833 nm (a) and at A = 800 nm (b). The pump 
pulse creates vibrational WPs in various electronic states. 
Significant ionization through the probe pulse occurs in (a) 
only if a WP is located around the outer, and in (b) only if a 
wave packet is near the inner turning point (FC window). 



collective degrees of freedom of the He droplets (phonons 
or ripplons). In our phenomenological approach we do 
not specify the microscopic dimer-droplet interaction. 
Instead, we use a well-established master equation, de- 
scribing damping for near-harmonic systems at effectively 
zero temperature fully quantum mechanically (see [22I]). 
Thus, the dynamics of the full vibrational density oper- 
ator p is determined from 



with earlier gas phase results Q. The weak overall de- 
crease of the oscillation amplitude in Fig. [2tb) is caused 
by the spreading of the initially well-localized WP in the 
anharmonic potential. This gas phase result is markedly 
different from the experimental HENDI spectrum (Fig. [21 
(a)), which is a clear sign of the effect of the He environ- 
ment. The measured signal shows an exponential decay 
of the initial amplitude which is well captured by our cal- 
culations (Fig. [21(c)). The persisting oscillations are also 
reproduced, yet with a smaller amplitude. 

Best agreement between experiment and our theory 
is obtained when damping, shift and desorption are in- 
cluded into the calculation with parameters 7 = 0.15/ps, 
A n = -50 cm" 1 and t d = 8 ps (Fig.H(c)). The depen- 
dence of all our findings on these precise values is rather 
smooth and only after leaving a ±5% interval significant 
deviations from the displayed figures occur. 

It is tempting to relate the damping rate 7 to mea- 
surements of the viscosity 77 in liquid He (see [16[). How- 
ever, such values for 77 are based on macroscopic theories 
which are not applicable here. Note that in [21| the au- 
thors investigate dissipative alkali dimer dynamics on a 
microscopic scale by studying collisions with 4 He atoms. 
Remarkably, the calculated friction coefficient for singlet 
systems (0.06/ps) turns out to be of the order of magni- 
tude of our 7. 

The temporal evolution and intensity of certain fre- 
quency components can be visualized by employing 
sliding window Fourier transforms (FT) on the time- 
dependent signal ("spectrogram" [§]). As discussed in 



Ref. 11( , in the experiment at 833 nm excitation wave 
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Here, H denotes the vibrational Hamiltonian and a,i,a\ 
are the annihilation/creation operators obtained from a 
harmonic approximation of the PES corresponding to 
electronic state i. We choose the parameters for damp- 
ing (rates ji = 7) and desorption (time td) to be state- 
independent. Note, however, a possible influence of the 
orientation of the dimer on damping dynamics as dis- 
cussed for LA2 in [2l|. We here assume no change of 
alignment during dynamics. 

Excitation at a wavelength of 833 nm almost exclu- 
sively follows the path sketched in Fig. Q] (a) and probes 
the WP propagation in the A 1 !]^ electronic state ficj . 
The probe pulse leads to an enhanced population of the 
final ionic state only if the WP is located at the outer 
turning point, where the intermediate 2n state opens a 
Franck-Condon (FC) window. The resulting oscillatory 
structure (Fig. (2)) in the ion rate reflects the dynamics of 
the WP in the A state. The simulation of the ion sig- 
nal for the gas phase is shown in Fig. [21b) and agrees 
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FIG. 2: Pump-probe signal for A = 833 nm. (a) Experimen- 
tal HENDI result (from Ref. [H]]); (b) Numerical simulation 
for the gas phase (see also 0|); (c) Calculated signal assum- 
ing He-induced damping, shift of PES and desorption of the 
molecule. 
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length a small increase of the frequency of the WP mo- 
tion during the first 10 ps is observed. In Fig. [3] this 
behavior (if less pronounced) is seen in the FT spectrum 
of our calculated signal, provided damping is included. 
The findings allow the following interpretation: during 
the first few picoseconds the contribution of the damped 
molecules is most important. As a consequence of the an- 
harmonic potential, they vibrate with a higher frequency, 
explaining the initial shift of the main frequency o;^ 33nm 
to slightly larger values (or lower vibrational quantum 
numbers). Later on, for times r > 10 ps, vibrationally 
damped molecules no longer contribute, as a closing of 
their FC window takes place. Hence, only molecules 
which desorb very early and thus do not suffer vibra- 
tional damping contribute to the ion yield at late delay 
times. Consequently, a weak oscillatory signal persists 
in the time domain and in the windowed FT the main 
frequency returns to the original gas phase value. 

At the shorter wave length A = 800 nm the excita- 
tion scheme sketched in Fig. [1] (b) applies. Here, a WP 
in the ground state, excited through stimulated Raman 
scattering (FJSRS), can be visualized through 3-photon 
absorption. In the gas phase, this process is superposed 
coherently with transitions from higher lying electronic 
states near their respective inner turning points, lead- 
ing to interferences in the ionic state. A more detailed 
analysis shows that in the gas phase the interference is 
indeed destructive: the contribution of the WP on the 
X surface is missing in t he g as phase (see Fig. 2] (b)) 
- see similar findings in [231 ]. Without any damping 
mechanism, contributions from the A state (component 
w 800nm ra 64cm" 1 ) and from the 2H state (component 
2w| n 0nm 86cm _1 ) are dominant, which is not observed 
in the He droplet experiment (Fig.0|a)). There, although 
w™™ is present during the first 5 ps, it fades out and 
only the component u>^ 0nm « 91 cm -1 (WP in X state) 
contributes at later delay times. We now include damp- 
ing for A = 800 nm with precisely the same parameters as 
for the considerations at A = 833 nm. The correspond- 
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FIG. 3: Calculated pump-probe signal for A = 833 nm. The 
dominant frequency is slightly shifted to larger values for short 
delay times (the black solid line indicates the the mean fre- 
quency between the dashed lines as a function of time). 
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FIG. 4: Spectrograms showing the WP dynamics at A = 
800 nm for the experimental HENDI result (a), for the gas- 
phase simulation (b), for the fully damped model (c), and for 
simulated data including a critical velocity (d). 



ing spectrogram (Fig. [He)) agrees remarkably well with 
experimental findings in many respects. The component 
at w^ 00nm fades out after several ps due to vibrational 
damping: the WP in the A state leaves the initial FC 
window at the inner turning point. 

More significant is the behavior of the component near 
w™™ 1 , ascribed to the ground state WP. It is clearly vis- 
ible in the experimental and theoretical (including damp- 
ing) spectrogram, yet absent in the gas phase. Damping 
of the WPs in the A and 2n state enhances ionization 
from the ground state by the probe pulse, since the map- 
ping of the X state dynamics to the ion no longer suffers 
from destructive interference with competing processes 
near the inner turning point on the A and 2n surfaces. 
Most remarkably, the results at 800 nm hint at a direct 
influence of superfluidity on these spectra. Molecules are 
assumed to be able to move unhindered through a su- 
perfluid as long as their velocity does not exceed the 
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FIG. 5: Dynamics in velocity domain of a WP excited in the 
electronically excited A-state and in the ground state X. The 
hatched area indicates velocities below the critical value of 60 
m/s. 



(Landau) critical velocity v c [B|. In a simple approach 
we determine from our calculations the average velocity 
(v) = (P)/2/x of the nuclear WP. Once it drops below the 
critical velocity (we chose the bulk value v c = 60 m/s), 
we set the friction coefficient to zero (see also 24J) and 
compare with the fully damped model. 

Fig. [5] illustrates the WP velocities in the ground (X) 
and excited state (A) . For the excitation at A = 833 nm, 
it is the dynamics on the A surface that dominates the 
observations. This vibration, however, is fast and the 
corresponding velocity oscillates between ±500 m/s. For 
most of the time, therefore, friction is present and the 
influence of v c is hardly noticeable. For an excitation at 
wave length A = 800 nm, however, the motion on the X 
surface is so slow (at most v « 15 m/s) that we neglect 
damping on that surface entirely (jx = 0). In the re- 



sulting spectrogram (Fig. EUd)), the component w|P 0nm 
is even more pronounced compared to the previous cal- 
culation (full damping). Most strikingly, the agreement 
with the experimental result is improved. 

In conclusion, the study of vibrational wave packets in 
dimers attached to He nanodroplets exhibits clear devia- 
tions from previous gas phase results. A model based on 
shifts of potential energy surfaces, desorption from the 
host and, most importantly, damping of the wave pack- 
ets is able to explain crucial features of the measurement. 
For excitation at A = 833 nm the decay of the signal and 
a small shift in the FT spectrum are reproduced. Damp- 
ing of excited state wave packets enhances the mapping 
of the ground state dynamics to the ion signal, as realized 
at A = 800 nm. 

Most remarkably, in our studies best agreement with 
experiment is obtained when we neglect damping for very 
slowly moving wave packets. We attribute this observa- 
tion to the role of the critical velocity in these experi- 
ments, which turns out to be larger than the velocity of 
the motion on the X surface. Thus, the influence of su- 
perfluidity is directly visible, opening the door for further 
studies: choosing heavier molecules for which the vibra- 



tional velocities are well below some 100 m/s, in combi- 
nation with a microscopic fully quantum calculation will 
enable us to investigate friction and frictionless motion 
on the atomic scale. 
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